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Abstract 



A model for the pionic components of nuclear wave functions is obtained 
from light front dynamical calculations of binding energies and densities. The 
pionic effects are small enough to be consistent with measured nuclear di- 
muon production data and with the nucleon sea. But the pion effects are 
large enough to predict substantial nuclear enhancement of the cross section 
for longitudinally polarized virtual photons for the kinematics accessible at 
Jefferson Laboratory, 
shadow / pion / reveal.tex 
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The belief that pions are the carrier of the nuclear force and hence the dominant agent 
responsible for the binding of nuclei has been a basic premise of nuclear physics since the 
time of Yukawa's original theory. This tenet was strengthened via the observation of very 
significant effects of pion exchange currents in a variety reactions involving electromagnetic 
probes of nuclei ]1J. Furthermore, chiral Lagrangians involving pions and nucleons are 
believed to be the correct low-energy representations of QCD 0. 

All of this basic physics was brought into question by lepton-nucleus deep inelastic scat- 
tering experiments and related measurements of the production of Drell-Yan fi + fi~ pairs in 
high energy proton-nucleus interactions. The lepton-nucleus deep inelastic scattering ob- 
served by the EMC experiment |3|.|4j] showed that there is a significant difference between 
the parton distributions of free nucleons and nucleons in a nucleus, including the substantial 
result that the nuclear structure function is supressed in the region of the Bjorken variable 
x ~ 0.5. This means that the valence quarks of bound nucleons carry less plus-momentum 
than those of free nucleons. But partons must carry the total momentum of the nucleus so 
the plus-momentum of non-nucleonic constituents must be enhanced. One useful postulate 
was that it is nuclear pions which carry a larger fraction of the plus-momentum in the nu- 
cleus than in free space |)]||. Such models [|j] explain the shift in the valence distribution 
and provide also an enhanced nuclear anti-quark distribution, observable in Drell-Yan ex- 
periments J7|. However, no such enhancement was observed ||, and no substantial pionic 
enhancement was found in (p,n) reactions 0. These failures to observe the influence of 
nuclear pions were termed a severe crisis for nuclear theory JlO . 



But the magnitude of the crisis depends on the size of the pionic effects expected from nu- 
clear theories. The earliest calculations [H] included enhancements via the effects of ttN <-> A 
transitions incorporated by summing the RPA series. Another calculation, involving signif- 
icant A effects, evaluated the expectation value of the pionic number density operator in 
a variational wave function [|ll|] and the result was characterized as infinite nuclear matter 
containing about 0.18 excess pions per nucleon. Both sets of calculations PJTT1] led to sig- 
nificant nuclear pion content and a corresponding enhancement of the nuclear anti-quark 
distribution, which was not observed experimentally [§]. However, it was recognized al- 
most immediately that one could reduce [I2| the predicted sea enhancement by a very small 
increase in the parameter g' used by || to represent the short-ranged nucleon-nucleon cor- 
relations. Thus it seems desirable to make further calculations of nuclear wave functions in 
a manner which includes pionic components along with nuclear saturation properties. 

Furthermore, it is natural to expect that calculations of high energy observables would 



benefit by using light front dynamics ||13|| . These considerations led us to use light front 
dynamics to compute nuclear wave functions WM. One approach was to use a new light- 
front one boson exchange potential to compute the saturation properties of infinite 
nuclear matter in a Bruckner theory calculation which obtained very reasonable values of 
the binding energy density, and compressibility. This calculation goes beyond mean field 
theory, as is necessary to obtain a non-zero pion content for nuclei with equal numbers of 
neutrons and protons. Another result is that the number of excess pions per nucleon is 0.05 
which seems to be in the range shown [12| as allowed by the Drell-Yan data. This smaller 
value of the pion excess arises mainly from the use of one-boson exchange, in which the A is 
absent ||16|| . Thus a reasonable nuclear theory exists which could coexist with the di-muon 



production data, but it is highly desirable to verify this set of dynamics by computing the 
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di-muon production cross section and by observing a non-zero signal of pionic effects. These 
are the aims of the present work. 

One possibility for verification arises from measurements |I7j of the ratio R, of scattering 



of virtual photons in a longitudinal (L) or transverse polarization state (T), R = o"l/o"t from 
nuclei. A large value of ol, and the corresponding violation of the Callan-Gross relation 



indicates the presence of nuclear bosons as fundamental constituents of nuclei |19|| . Indeed, 
the HERMES collaboration found a factor of five enhancement of R, but this was for values 
of x and Q 2 at which the pionic effects are too small to be relevant |jl9 . 



Here we study the possible efffects of nuclear pions on <tl as a function of x and Q 2 and 
therefore to determine if such effects are detectable. We shall also compute those effects 
for the nucleon target to predict its the pionic content. These effects are also relevant to 
experiments aimed at determining the pion elastic form factor [00] . 

The start is to confront the deep inelastic scattering data on the nucleon and the nuclear 
di-muon production data which place severe constraints on the pionic content of the nu- 
cleon [21] and nuclei, and ensure coexistence of the dynamical model with those data. The 



necessary phenomenology has been described often, see the reviews [|J. The main point is 
that the quark and anti-quark distributions can be given as convolution of the q, q distri- 
butions in a given nuclear hadronic constituent with the light cone distribution functions: 
fir/A(y), fn/N(y), the probability to find an excess pion in the nucleus (A) or nucleon (N), 
with a plus-momentum given by ym^. These quantities are obtained from the square of the 
ground state nuclear of nucleon wave function, computed using light front dynamics, as: 

U/n(v) =m N J d 2 k L (N\n n (k)\N) (1) 

U/a(v) = jrn N J d 2 k ± (A\ n7r (k)\A} - U /N (y) (2) 

where n n (k = (kj_,A; + )) is the pion number operator. The function f n /N{y) is constrained 



by the data on the nucleon sea which restrict the u and d distributions to be similar [2T 
The important input is the 7rN form factor, which we take to be of the Cloudy Bag Model 
form P21 with a bag radius, 0.9 fm, large enough to be consistent with the constraints. 



The quantity f n /A(y) controls the nuclear pionic effects. We would like to obtain this 
directly from our latest light front calculation [15 1, of nuclear matter properties. This calcu- 



lation provides the integral quantity for the excess number of pions as J °° dy f n /A(y) = 0.05, 
but not the function f w /A itself. However in all models this function is expected to be van- 
ish for very large and very small values of y and smoothly rise to a maximum, at about 
y = 0.2, 0.3. Thus the overall strength of the effects in Drell-Yan scattering and o"l we 
wish to examine here are not expected to be very sensitive to details or have much model 
dependence as long as the overall normalization is constrained. Therefore we revive our old 
calculation [JT3], but now using parameters which yield 0.05 for the number of excess pions 
in nuclear matter. The simplest way to achieve this value is to vary the value of g'. The 
relation g' = 0.795 gives the desired value for nuclear matter and also that excess nuclear 
pions carry a fraction 0.016 of the nuclear plus momentum. We account for the difference 
between Fe and nuclear matter by reducing the Fermi momentum from 1.37 fm -1 to 1.3 



fm [23]. This reduces the excess pion number to 0.04 and the momentum fraction to 0.135. 



The results for the pionic distribution functions are shown in Fig. p]. 
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FIGURES 
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FIG. 1. Excess pion distribution^functions for nuclei (solid) and for the nucleon (dashed). The 
results for different nuclei 12 C, 56 Fe, and 197 Au are obtained from the Fermi gas model with values 
of k f = 1.2, 1.3 and 1.34 fin" 1 [22]. 

We see that f n /N,A differ from zero for y ~ 0.2 — 0.5, but f n /A vanishes rapidly as y 
approaches 0. This means that pionic effects are irrelevant to understanding the HERMES 
effect, which occurs for xbj ~ 0.01. 
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FIG. 2. Ratios of the Drell-Yarr di-muon ratio per nucleon, for positive The data 

are from Ref. [8]. 

Using fv/A one can compute the contribution to the di-muon production cross section in 
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which a quark from the incident proton is annihilated by a nuclear anti-quark (from either 
a nucleon or a pion) to form a virtual photon which emerges as a di-muon-/i + /i~ pair. The 
pionic contributions must be added to the contributions of the nucleons. The nucleonic 



effects are computed using the formalism of p4j and [25|. The key parameter is the average 



separation energy, which is take as 56 MeV, as obtained from an extrapolation of the nuclear 
matter value of 71 MeV. Such parameters allow a reasonable qualitative representation of 
the nuclear deep inelastic scattering data P, pB] , pB |. 



With the present model, it is the nucleonic effects which control the computed Drell-Yan 
results. These are shown in Fig. Q for values of the target plus momentum fraction, x t 
large enough so that the effects of shadowing are absent. The compuations are made using 
the experimental acceptance f|. We see that it is possible to reproduce the experimental 
data with a theory which is based on detailed nuclear dynamics such as light front Bruckner 



theory |15 



The question remaining is, "Can there be a signal to observe these pions?" Thus we turn 
to the role of pionic effects on for nuclear and nucleon targets. The pionic contribution 
to the hadronic tensor is given by: 

FW"' = / A (P| J£(£) ^(0)|P), (3) 

where \P) represents the nucleon or nuclear target ground state of total momentum P, and 
the momentum of the virtual photon is q = [y, 0^, —y/u 2 + Q 2 ). The current operator which 
accounts for the effects arising from the nuclear pions is 

J£ = i ((f)*d^(f) -<f)dy*), (4) 

with as the complex pion field operator. The effects of the pion charge form factor F n (Q 2 ) 
are included below. 

We are concerned with <tl, and use the standard general formula p7| , p^| ; 

° L ~ (Q 2 + V 2 )^ 2 W ■ W 

We now evaluate 5 n W m by using J° in Eq. (Q). Then insert a set of states which asymp- 
totically contain a nucleus and a free pion of momentum p^. The sum over nuclear states 
can be performed using closure if one uses light front variables (q ± = q° ± q 3 ) to describe 
the momentum of the photon and of the nuclear pion (kj_, k + ). For large enough values of 
v the interactions of the struck pion with the residual nucleus may be ignored, the energy 
of the outgoing pion may be taken as z/, and the effects of non-zero values of nuclear excita- 
tions energies (proportional to q + /q~) appearing in the energy-momentum conserving delta 
function inherent in Eq. @ may be ignored. Then the result of a straightforward evaluation 
leads to the result: 

S-W™(N, A) = ~ U/{ "f {0 F*(Q\ (6) 



2Q 2 3 M 



N 



where F n (Q 2 ) is taken as a monopole form consistent with the observed mean square radius, 
and the factor of 2/3 accounts for the fact that only charged pions enter in electron scattering 
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(with N=Z). Note that energy momentum conservation gives the relation k 
and £ is the Nachtmann variable. 
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FIG. 3. Pionic contribution to the photon longitudinal cross section on the proton, as a function 
of Q 2 and x. 

The pionic contribution to the longitudinal cross section for a proton (p=N) target, 
a£(p). is obtained from Eqs. @ and (|[), and is displayed in Fig @. Cross sections of this 
size are routinely observable. The decrease with increase of Q 2 is caused by the pion form 
factor. The value of crj(p) for x = 0.2, Q 2 = 0.7 GeV 2 corresponds to a contribution to R 
of the proton of about 0.04, small enough to avoid any contradiction with existing data. 

We now turn to the case of nuclear targets, The standard general formula (H) allows 
us to relate 5 ,r VF 00 to the corresponding contribution (per nucleon) d^a^ to the longitudi- 
nal cross section, so that the nuclear longitudinal cross section (per nucleon) is given by 
&l(A) = ol(D) + 8* a i,. It is desirable to present a ratio a L {A) j a L {D\ or equivalently 
(Jl(A) / ctt(D)R D) which allows the use of parameterizations for F 2 (D) and R(D). Thus we 
obtain the result: 



1 + 



5 n W 



oo 



(Q 2 + v 2 )vAF 2 D R D 



;i + Rd) 



(7) 



which may be evaluated using Eq 



to be: 



2 v 2 
a L (D) - 1 + X 3 m) (Q 2 + u 2 \ 



F 2 D R D 



A + Rd) 



(8) 



where as usual x = Q 2 ' jlM^v. Note that only the pionic effects are included here. The 
effects of nuclear vector and scalar mesons, so important in understanding the HERMES 



effect [|19|], are not included because these appear at much lower values of x than we shall 
consider. 
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The results of our calculation of <tl on 56 Fe are shown in Fig. f| which shows the difference 
between the left hand side of Eq. @ and unity. Results for other nuclei can be obtained by 
using the different pion distribution functions of Fig. [I]. 
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FIG. 4. Enhancement of longitudinal cross sections, as a function of Q 2 and x. 

Large effects are predicted for a quantity which is readily measurable at Jefferson Lab- 
oratory. The largest effects obtained with the lowest value of Q 2 = 0.3 GeV 2 , correspond 
to very low values of v ~ 0.4 GeV, where there may be substantial corrections to the light- 
front closure approximation because q + is not vanishingly small. While there seems to be 
no reason why the effects of such corrections would cancel the effects shown here, the inter- 
pretation of the data would be simplified by observing effects at larger values of Q 2 . In this 
case one obtains smaller values of g + , and larger values of z/, but the outgoing pion energies 
are in the vicinity of the 7r-nucleon resonance region, so that final state interaction effects 
could change the shape of the curve. However, since the variable <jl involves an inclusive 
measurement the overall strength as represented by, e.g. the area under the curve, should 
not be affected very much. 

The substantial 15-30% effects, for values of x ~ 0.2 — 0.3, present an opportunity to 
observe a significant nuclear enhancement of a^. This is an excellent opportunity to unravel 
a significant long-standing mystery involving the absence of nuclear pionic effects. 
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